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Abstract 

The Busse-Heikes dynamical model is described in terms of relaxational and nonre- 
laxational dynamics. Within this dynamical picture a diverging alternating period 
is calculated in a reduced dynamics given by a time-dependent Hamiltonian with 
decreasing energy. A mean period is calculated which results from noise stabiliza- 
tion of a mean energy. The consideration of spatial-dependent amplitudes leads 
to vertex formation. The competition of front motion around the vertices and the 
Kiippers-Lortz instability in determining an alternating period is discussed. 



1 Introduction 



One of the most extensively studied systems, in the field of pattern formation 
in nonequilibrium systems, is Rayleigh-Benard thermal convection. In many 
geophysical and astrophysical systems, thermally induced convection is com- 
bined with Coriolis forces induced by rotation. Therefore, Rayleigh-Benard 
convection in fluid layers rotating around a vertical axis is a hydrodynami- 
cal system of significant importance. Specially interesting is a spatio-temporal 
regime that takes place above a critical rotation angular velocity. The system 
breaks up into a persistent dynamical state such that set of parallel convec- 
tion rolls are seen to change orientation with a characteristic period. This 
phenomenon is known as the Kiippers-Lortz instability [1]. This instability 
can be described as follows: for an angular rotation speed Q greater than 
some critical value Q c , convective rolls lose their stability with respect to rolls 
inclined at an angle of about 60° in the sense of rotation. The new rolls undergo 
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the same instability, so that there is no stable steady-state pattern. As a re- 
sult spatially disordered patterns arise already arbitrarily close to the onset of 
convection. Experimental characterization of this regime of spatio-temporal 
chaos has been reported in [2]. Most experiments have been performed for 
small Prandtl numbers and have been theoretically described in the realm of 
Swift-Hohenberg models with [3] and without [4] the inclusion of mean flow 
effects. On the other hand, large Prandtl numbers lead to more rigid convec- 
tion rolls. In this situation mean flow coupling can be neglected and, in the 
limit of infinite Prandtl numbers, three-mode models have been shown [5,6] 
to exhibit the same qualitative features as more sophisticate Swift-Hohenberg 
models that take into account the full range of possible roll orientations. 

We consider in this paper a three-mode model proposed by Busse and Heikes 
[7] to study the Kiippers-Lortz instability. Each mode represents the ampli- 
tude of a set of parallel rolls with an orientation of 60° to each other. This 
model contains an attracting heteroclinic cycle connecting three fixed points 
corresponding to the three different roll solutions. The model predicts success- 
fully the existence of a region in parameter space in which the roll solution is 
unstable, but fails to reproduce the experimental observation of an approxi- 
mately constant period between roll alternation. Whereas Busse and Heikes 
speculated that such a constant period would be obtained by the addition of 
noise, a conclusion confirmed by Stone and Holmes [8], no systematic study 
of the relation of the period to the system parameters has been performed 
so far. Another explanation for period stabilization has been given by Cross 
and Tu [5] who have performed numerical investigations of an extension of the 
Busse-Heikes equations, where a spatial variation of the amplitudes has been 
introduced. In this paper, we study in detail these two proposed mechanisms 
for period stabilization in the Busse-Heikes model: (i) addition of noise and 
(ii) the consideration of spatial-dependent terms. 

The paper is organized as follows: in section 2 we present a description of 
the Busse-Heikes model and give a clear physical explanation of the period 
divergence. We describe the dynamics in terms of a relaxational and a nonre- 
laxational part. The alternating period is calculated for the latter part which 
is associated with a slowly varying time-dependent Hamiltonian. In section 3 
we consider the same model with the inclusion of additive white noise terms 
and we calculate the mean period stabilized by noise in terms of the previous 
dynamical picture. Sections 2 and 3 discuss ordinary differential equations for 
the amplitudes of the three modes. In section 4 we consider the more physi- 
cally appropriate situation of spatial-dependent amplitudes in a d = 2 model 
and study the influence on the dynamics of isotropic and anisotropic spatial- 
dependent terms. We describe the formation of vertices and how the period of 
roll alternation is determined by the competition of front motion around the 
vertices and the Kiippers-Lortz instability. 
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2 Busse-Heikes Model 



Based on the fact that, in a first approximation, only three directions are 
relevant to this problem, Busse and Heikes [7] proposed a dynamical model to 
study the Kiippers-Lortz instability. The vertical component of the velocity 
field ^(r, t) is written as: 

3 

*(r, t) = M T > t)e iqo ^- r + c.c. (1) 

3=1 

("c.c." denotes complex conjugate). The vectors ej are unit vectors in direc- 
tions j = 1,2,3 which form an angle of 60° between them, and go is the 
selected wavenumber of the convection pattern. In this model the (complex) 
amplitudes of the three rotating modes, Ai, A 2 , A 3 , are independent of space 
and follow the evolution equations [7]: 



M = A 1 [u- \A,\ 2 - (1 + n + 5)\A 2 \ 2 - (1 + n - 8)\A 3 \\ 

A 2 = A 2 \v - \A 2 \ 2 - (1 + /i + b)\A 3 \ 2 - (1 + fj, - 5)^1% (2) 

A 3 = A 3 [v - \A 3 \ 2 - (1 + fj, + 5)|Ax| 2 - (1 + n - 5)\A 2 \ 2 }. 

The parameter v is proportional to the difference between the Rayleigh num- 
ber and the the critical Rayleigh number for convection. We will consider 
exclusively in this paper the case of well-developed convection for which the 
parameter v can be rescaled to 1, i.e. v — 1 henceforth. The exact relation of 
ji and 5 to the fluid properties has been given in [1]. We mention here that 
fi is a parameter related to the temperature gradient and the Taylor number 
(proportional to the rotation speed ft) in such a way that it takes a nonzero 
value in the case of no rotation, Q = 0, whereas 5 is related to the Taylor 
number in such a way that Q = implies 5 = 0. We will consider only fl > 0, 
or 5 > 0; the case f2 < (5 < 0) follows by a simple change of the coordinate 
system. Although the dynamical equations are defined for all values of the 
parameters, only the case \x > is physically relevant. Writing Aj = ^Rje l8j 
we obtain equations for the modulus square of the amplitudes Rf 



fli = 2i2i[l - i?i - (1 + n + 5)R 2 - (1 + n - S)R 3 ], 

R 2 = 2R 2 [l-R 2 -(l+/2 + 5)R 3 - (1 + - 5)R 1 ], (3) 

R 3 = 2R 3 [l-R 3 -(l+fi + 5)R 1 -{I S)R 2 ], 

and for the phases 9f 
01 = 0, 
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#2 = 0, 
#3 = 0. 



(4) 



It follows that the phases are simply arbitrary constants fixing the location of 
the rolls. A solution of the form ^(r) = ^jR~e^ qi,e ^' r+e ^ + c.c. represents a set 
of rolls of wavelength 27r/g , oriented in a direction perpendicular to the vector 
&j. Hence, in this model one can simply consider the equations for the real 
variables Rj instead of the equations for the complex variables Aj. A similar 
set of equations has been proposed to study population competition dynamics. 
For a single biological species, the Verhulst or logistic model assumes that its 
population N(t) satisfies the evolution equation: 

^L = rN(l-\N), (5) 



where r is the reproductive growth rate and A is a coefficient denoting com- 
petition amongst the members of the species. If three species are competing 
together, it is adequate in some occasions to model this competition by intro- 
ducing a Gause-Lotka-Volterra type of equations [9,10]: 
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which are the same that the Busse-Heikes equations (3) for the modulus square 
of the amplitudes Rj with the identifications: r — 2, X — 1, a — l + f/, + 5, 
j3 = 1 + /i — 5. These equations are the basis of May and Leonard analysis 
[9]. We also mention the work of Soward [11] which is concerned with the 
study of the nature of the bifurcations mainly, but not limited to, close to 
the convective instability for small v, in a slightly more general model that 
includes also quadratic nonlinearities in the equations. In the remaining of the 
section we will analyze some of the properties of the solutions of the Busse- 
Heikes equations (2). Although our analysis essentially reobtains the results 
of May and Leonard, we find it convenient to give it in some detail because, 
besides obtaining some further analytical expressions for the time variation of 
the amplitudes, we are able in some cases of rewriting the dynamics in terms 
of a Lyapunov potential. The existence of this Lyapunov potential allows us 
to interpret the asymptotic dynamics for \i = as a residual (conservative) 
Hamiltonian dynamics. For \i > we will use an adiabatic approximation 
with a time-dependent Hamiltonian. This interpretation will turn out to be 
very useful in the case that noise terms are added to the dynamical equations, 
because the found Lyapunov potential governs approximately the stationary 
probability distribution. 
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We first look for stationary solutions of the Busse-Heikes equations (2). The 
fixed point solutions are the following: 



(a) The null solution: R ± = R 2 = R 3 = 0. 

(b) Roll solutions. There are three families of these solutions, each character- 
ized by a unique nonvanishing amplitude, for instance: R 2 , R3) = (1, 0, 0) 
is a roll solution with rolls perpendicular to the ei direction, and so on. 

(c) Hexagon solutions. The three amplitudes are equal and different from 0, 
namely R\ = R 2 = R3 = They only exist for ji > —3/2. 

(d) Rhombus solutions. There are three families of these solutions, in which 
two amplitudes are different from and the third amplitude vanishes. For 
instance: (R ± , R 2 , R 3 ) = ( -^M_ T > , 0). They only exist for fi > 5, or 

-1 - VTTP < /2 < -5. 

The stability of the previous solutions can be studied by means of a linear 
stability analysis. The result is summarized in Fig. 1. For /i < —3/2 there are 
no stable solutions and the amplitudes grow without limit. The rhombus and 
null solutions are never stable. The hexagon solutions are stable for — 3/2 < 
/i < 0. The roll solutions are stable for /i > 5. For < /i < 5 there are no 
stable solutions, but the amplitudes remain bounded. This instability can be 
described as follows: consider the unstable roll solution R 2 , R3) = (1, 0, 0). 
The amplitude of the A 2 mode starts growing and that of A\ decreasing in 
order to reach the roll solution (0, 1, 0). However, this new roll solution is also 
unstable, and before it can be reached, the dynamical system starts evolving 
towards the roll solution (0,0, 1), which is unstable and evolves towards the 
solution (1, 0, 0) which is unstable, and so on. Schematically, we can represent 
the situation as follows: 



This is the Kiippers-Lortz instability that shows up in the rotation of the 
convective rolls. The Kiippers-Lortz unstable region is characterized by the 
presence of three unstable fixed points, and a heteroclinic cycle connecting 
them. 

The novelty of our treatment consists in writing the Busse-Heikes equations 
of motion in the form: 



(1, 0, 0) - (0, 1, 0) - (0, 0, 1) - (1, 0, 0) - (0, 1, 0) . . . 



(7) 




dV 
dA* 



+ Svj, 



J = 1,2,3, 



(8) 



with the potential function: 
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V(A 1 ,A 2 ,A 3 ) 
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(9) 



= -(R 1 +R 2 + R 3 ) + 1 (i? 2 + ^2 + ^ + (1 + ^ (fl^ + # 2 i? 3 + ^ 3jRi ) ? 



«1=Ai( 


-\A 2 \ 2 + \A 3 \ 2 ) 


= A 1 (-R 2 + R 3 ), 




v 2 = A 2 ( 


-\A 3 \ 2 + \A,\ 2 ) 


= A 2 (-R 3 + R 1 ), 


(10) 


V3 = A 3 { 


-\A 1 \ 2 + \A 2 \ 2 ) 


= A 3 (-R 1 + R 2 ). 





The first term in the right-hand side of (8) describes relaxation in the potential 
V(A 1 , A 2 , A 3 ). In the case 5 = 0, hence, the dynamics is described simply as 
the relaxation, along the gradient lines of the potential V, in order to reach 
a minimum of V. In the case 5 > there is another contribution to the 
dynamics. Its effect can be analyzed partly by looking at the time evolution 
of the potential: 



~dt ~ dA~ Hi + C ' C ' = 



dV 



OA, 



+ 5 



j=l OA 3 



■ (11) 



Therefore, when the so-called orthogonality condition is satisfied: 



^dA V] + c - c - 



o, 



(12) 



the function V decreases along the dynamical trajectories and it becomes a 
Lyapunov potential [12] if it is bounded from below (which is the case for 
pi > —3/2). Using equations (9) and (10), (11) can be rewritten as: 

- 2^(|A| 2 - \A 2 \ 2 )(\A 2 \ 2 - \A 3 \ 2 )(\A 3 \ 2 - |A| 2 ),(13) 



dV_ 
~dt 



-2E 



dV 



dA j 



so that the orthogonality condition is seen to be satisfied for /j,S = 0. In 
the case 5 = 0, the system is purely relaxational in the potential V and 
the corresponding stability diagram can be obtained also by looking at the 
minima of V. For the null solution, the potential takes the value V = 0; 
for the rhombus, V = —1/(2 + //); for the roll solution, V = —1/2; and, 
finally, for the hexagon solution, V = —3/(6 + 4//). The study of the potential 
(for 5 = 0) shows that the null and rhombus solutions correspond always to 
maxima of the potential and are, therefore, unstable everywhere. It turns out 
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that the rolls (hexagons) are maxima (minima) of the potential for fi < 
and minima (maxima) for /i > 0. Also, the potential for the roll solution is 
smaller that the potential for the other solutions whenever /i > 0, indicating 
that the rolls are the most stable (and indeed the only stable ones) solutions 
in this case. Unfortunately, this simple criterion does not have an equivalent 
in the nonrelaxational case, 5 > 0, for which one has to perform the full linear 
stability analysis. 



2.1 The case /i = 

According to the result (13), the function V(A 1 , A 2 , A 3 ) is a Lyapunov poten- 
tial whenever fi5 = 0. As discussed in the previous section, the case 5 = 
implies a relaxational gradient dynamics in which all variables tend to fixed 
values. In the case // = 0, 5 > 0, the dynamics is nonrelaxational poten- 
tial [13-15] and, whereas the dynamics still leads to the surface of minima of 
the Lyapunov function, there is a residual motion in this surface for which 
dV/dt = 0. In other words: the relaxational terms in the dynamics make the 
system evolve towards the degenerate minimum of the potential (which for 
\i = occurs at R± + R 2 + R 3 = 1). The residual motion is governed by 
the nonrelaxational part which is proportional to 5 and this residual motion 
disappears for 5 = 0, the relaxational gradient case. 

According to this reduction of the dynamics as a residual motion in the surface 
of minima of the potential V, strictly valid only for /i — 0, it turns out that 
it is possible to solve essentially the equations of motion. By "essentially" we 
mean that after a transient time in which the system is driven to the minima 
of V, the residual motion is a conservative one in which it is possible to define 
a Hamiltonian-like function that allows one to find explicit expressions for the 
time variation of the dynamical variables. Let us define the variable 



It is straightforward to show that, for arbitrary jj, and 5, X satisfies the evo- 
lution equation: 



X(t) = R 1 + R 2 + R 3 . 



(14) 



X = 2X(1 - X) - AfiY, 



(15) 



where: 



Y(t) = R 1 R 2 + R 2 R 3 + R 3 R 1 . 



(16) 
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In the case fj, — the equation for X(t) is a closed equation whose solution is 

X(t) = t 1 . (17) 



Here Xo = X(t = 0). From this expression it turns out that lim t ^ 00 X(t) = 1 
independently of the initial condition. In practice, and due to the exponential 
decay towards 1 of the above expression, after a transient time of order 1, X(t) 
already takes its asymptotic value X(t) = 1. Therefore, we can substitute 
Ri(t), say, by 1 — R2(t) — Rs(t) to obtain evolution equations for R2{t) and 
Rs(t). In this way, the original 3- variable problem, Eqs. (3), is reduced to a 
residual dynamics in a 2- variable subspace: 



R 2 = 25R 2 (1-R 2 -2R 3 ), (18) 

R 3 = -25R 3 (1-2R 2 -R 3 ). (19) 

These are Hamilton's equations: 

A = 2*J* (20) 

corresponding to the Hamiltonian: 

H(R 2 , Rs) = R 2 R 3 (l -R2- Rs)- (22) 



As a consequence, in the asymptotic dynamics for which the Hamiltonian 
description is valid, H(t) is a constant of motion, H = E, which will be called 
the "energy" . The Hamiltonian dynamics is valid only after a transient time, 
but the value of E depends only on initial conditions at t — 0. The dependence 
of E on the initial conditions can be found by introducing the variable 7i: 

1~L = R\R 2 R 3 (23) 



which, in the asymptotic limit (t — > oo) is equivalent to 7i. It is easy to 
show that, for arbitrary values of /i and 5, ji satisfies the following evolution 
equation: 

7i- 1 ^ = 6-(6 + 4/i)X (24) 
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(one can reduce the original dynamical problem to variables {X, Y, 7i} but 
the equation for Y turns out to be too complicated, see [11]). If we substitute 
the solution for X(t) valid in the case /x = we obtain: 

(25) 



H(t)=H (1 - X )e' 2t + X { 



with 7i = H(t = 0). The asymptotic value for H is 

E = lim Hit) = lim H(t) = ^ = ^(0)^(0)^(0) ( } 

t-oo w *->«, w X 3 (i2x(0) + i? 2 (0) + i? 3 (0)) 3 V ; 



Again, this asymptotic value is reached after a transient time of order 1. This 
expression suggests to define the time- dependent variable: 



whose evolution equation (again, for arbitrary /x, 5) is: 

^ = -4/i (x - 3-D E = -W(t)E. (28) 



Therefore, in the case /x = 0, E(t) = E is a constant of motion that co- 
incides, in the asymptotic limit when X — 1, with the numerical value of 
the Hamiltonian 7i. According to their definition, E(t) is a bounded function 
< E(t) < 1/27 and /(*) > for Rj > 0, j = 1, 2, 3. 

The problem in the case /x = can now be given an explicit solution. After a 
transient time (or order 1), the motion occurs on the plane Ri + R 2 + -R3 = 1, 
see Fig. 2. The motion is periodic because it corresponds to a Hamiltonian 
orbit with a fixed energy. The exact shape of the trajectory depends on the 
value of the energy E which, in turn, depends on initial conditions. More 
interestingly, the period of the orbit can also be computed. For this, we solve 
the evolution equation (again asymptotically) for, say, R 3 . By elimination of 
R 2 by setting H = E in Eq. (22): 

Ri = \ (l - Rs ± ^(1-R 3 ) 2 -4E/R^j , (29) 

we obtain a closed equation for R 3 : 

Rs = ±25^/R 2 3 (l-R 3 y-AER 3 . (30) 
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Let b and c be the return points, i.e. the solutions of 

R 3 (l - R 3 ) 2 - AE = 0, 



(31) 



lying in the interval (0,1) [16]. The three roots, a,b,c, of the above third- 
degree equation are real and two of them (the return points b, c) lie in the 
interval (0, 1). The explicit expression for the roots is: 



a — - 



c— - 



1 + cos - 

o 



1 + COS 
1 + COS 



9-2n 
3 

9 + 2n 



(32) 
(33) 
(34) 



where 



9 = arccos (54E - 1). 



(35) 



Integration of (30) yields the equation of motion for i?3(t): 

Mt) t 
r dx =26 r (3g) 

•> Jx(x — a)(x — b)(x — r\ J 



where we have chosen the initial time to to correspond to the minimum value 
when Rs(t) = c. The integral in the left hand side can be expressed in terms 
of the Jacobi elliptic function [17] sn[x|g], to yield: 

R 3 (t) = bc_ 

b + (c - b)sn*[6y/b(a-c){t - t )\q] 



where 

a(b — c) 
^ b(a — c) 



(38) 



The period of the orbit T can be expressed in terms of the complete elliptic 
function of the first kind K(q): 



10 



and i?3(t) can be written as: 



be 



(40) 



6+(c-6)sn 2 msl(t-t )\qY 



Finally, the evolution equations for the other variables are: 



R 1 (t) = R 3 (t-T/3), 
R 2 (t) = R 3 (t-2T/3). 



(41) 
(42) 



Summarizing, the behavior of the dynamical system in the case \i = can 
be described as follows: after a transient time (or order 1) the three variables 
Ri, i?2, -R3 vary periodically in time on the plane R\ + R2 + R3 = 1. When 
Ri decreases, R 2 increases, etc. The period of the orbit depends only on the 
initial conditions through a constant of motion E. The explicit expression for 
the period, Eq. (39), shows that the period diverges logarithmically when E 
tends to zero, namely 



and the amplitude of the oscillations A = b — c depends also on the constant 
E. When E tends to the amplitude approaches 1: 



All these relations have been confirmed by a numerical integration of the 
Busse-Heikes equations. In Fig. 3 we plot the time evolution of the amplitudes 
in the case fj, — 0, 5 — 1.3. In this figure we can observe that, after an initial 
transient time, there is a periodic motion (characteristic of the Kiippers-Lortz 
instability) well described by the previous analytical expressions. 

2.2 The case /i > 

Once we have understood the case \x = 0, we now turn to /i > 0. In this case, 
the function V is no longer a Lyapunov potential and we can not reduce the 
motion to a Hamiltonian one on the surface of minima of V. However, since 
the main features of the Kiippers-Lortz dynamics are already present in the 
case fi = one would like to perform some kind of perturbative analysis valid 
for small /i in order to characterize the Kiippers-Lortz instability. We exploit 
these ideas in order to develop some heuristic arguments that will allow us to 
make some quantitative predictions about the evolution of the system. 




(43) 



A = (1-2£ 1/2 ) x (1 + 0(E)). 



(44) 
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According to Eq. (28), one can infer that E(t) decreases with time in a char- 
acteristic time scale of order yT l . If fi is small, E(t) decreases very slowly and 
we can extend the picture of the previous section by using an adiabatic ap- 
proximation. We assume, then, that the evolution for /i > can be described 
by a Hamiltonian dynamics with an energy that slowly decreases with time. 
Hence, in reducing the energy, the system evolves by spiraling from a periodic 
orbit to another (similarly to a damped harmonic oscillator). Assuming this 
picture of a time-dependent energy E(t), the main features of the case \i = 
can now be extended. This model has several predictions: 

• After a transient time of order 1, the motion occurs near the plane i?i + 
R<i + R3 — 1. This is checked in the simulations as we can see in Fig. 4 where 
we plot the time evolution of the three amplitudes as well as their sum, in the 
case 5 = 1.3, /i = 0.1. 

• The period of the orbits is now a function of time. Since the energy de- 
creases towards zero, it follows from Eq. (43) that the period diverges with 
time. Moreover, it is possible to give an approximate expression for the time 
dependence of the period. By integration of equation (28), we obtain: 

E(t) = £(t )e~ 4/ ^o mdt ' » E(t )e- 4 ^- t0 \ (45) 

where we have approximated fit) by its asymptotic value f(t) = 1. Once we 
have the time evolution of the energy, we can compute the time dependence 
of the period by using T(t) = T(E(t)) as given by (39). For late times, the 
energy is small and the asymptotic result (43) leads to: 

T(t) = To + *£t. (46) 

This shows that the period increases linearly with time, in agreement with the 
results of [9] in which the residence period was shown to behave also linearly 
with time (although with a different prefactor). In order to check this relation, 
we have performed a numerical integration of Eqs. (2) and computed the pe- 
riod T, defined as the time it takes for a given amplitude to cross a reference 
level (taken arbitrarily as Rj = 0.5), as a function of time. The results for 
5 = {1.3, 3} and fi = {0.1, 0.01}, plotted in Fig. 5, show that there is a per- 
fect agreement between the theoretical expression and the numerical results. 

• The amplitude of the oscillations, as given by the return points A(t) = 
b{t) — c(t) is now a function of time. Using expression (44) with an energy 
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that decreases with time as in Eq. (45) we obtain that the amplitude of the 
oscillations increases with time, see Fig. 4, and that it approaches 1 in a time 
or order t ~ ji^ 1 . More specifically, we have: 

1 - A(t) = (1 - Ao)e" 2 ^. (47) 

In summary, for the case \i > 0, the period of the orbits, which is a function of 
the energy, increases linearly with time and the amplitude of the oscillations 
approaches 1. We characterize in this way the increase of the period between 
successive alternation of the dominating modes, see Fig. 4, as an effect of 
the Hamiltonian dynamics with a slowly decreasing energy. This prediction 
of the Busse-Heikes model for the Kiippers-Lortz instability is unphysical, 
since the experimental results do not show such an increase of the period. 
Busse and Heikes were fully aware of this problem and suggested that noise 
terms ("small amplitude disturbances"), that are present at all times, prevent 
the amplitudes from decaying to arbitrary small levels and a motion which is 
essentially periodic but with a fluctuating period is established. In the next 
section we study the effect of noise in the dynamical equations. 



3 Busse-Heikes model in the presence of noise 

In order to account for the effect of the fluctuations, we modify the Busse- 
Heikes equations by the inclusion of noise terms: 

A, = A 1 [l- \A,\ 2 - (1 + + 5)\A 2 \ 2 - (1 + - b)\A 3 \ 2 \ + 6(0, 

A 2 = A 2 [1- \A 2 \ 2 - (1 + fj, + 5)\A 3 \ 2 - (1 + fj, - 5)|Ai| 2 ] + 6(t), (48) 

A 3 = A 3 [l - \A 3 \ 2 - (1 + n + 5)|Ax| 2 - (1 + n - 5)\A 2 \ 2 ] + &(<)■ 

We take the simplest case in which the ^(t) are, complex, white-noise pro- 
cesses [18] with correlations: 

Ut)e j (t')) = 2e5(t-t')5 ij . (49) 

As mentioned before, and in the case of parameter values lying inside the 
Kiippers-Lortz instability region, noise prevents the system from spending an 
increasing amount of time near any of the (unstable) fixed points. The mecha- 
nism for this is that fluctuations are amplified when the trajectory comes close 
to one of the (unstable) fixed points of the dynamics and the trajectory is then 
repelled towards another fixed point [8]. Hence, a fluctuating, but periodic on 
average, trajectory is sustained by noise. Within the general picture developed 
in the previous section, the main role of noise for ji > is that of preventing 
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E(t) from decaying to zero. This can be understood in the following qualita- 
tive terms: when noise is absent, the dynamics brings the system to the surface 
of minima of V, where the dissipative terms act by decreasing the energy in 
a time scale of order /i" 1 , see Eq. (45). The inclusion of noise has the effect 
of counteracting this energy decrease that occurs in the surface of minima of 
V. As a consequence, E(t) no longer decays to zero but it stabilizes around a 
mean value (E) . By stabilizing the orbit around that one corresponding to the 
mean value (E), fluctuations in the residual motion stabilize the mean period 
to a finite value. In order to check this picture, we have performed numerical 
simulations of Eqs. (48) for small noise amplitude e, using a stochastic Runge- 
Kutta algorithm [15]. The numerical simulations, see Fig. 6, show indeed that 
the trajectories have a well defined average period (T). 

From a more quantitative point of view, and according to the previous pic- 
ture, we can compute the mean period (T), which in the purely Hamiltonian 
case was a function of E, see Eq. (39), by using the same function applied 
to the mean value of E, i.e. (T) = T((E)). This relation has been checked 
in the numerical simulations. In Fig. 7 we plot the mean period (T) versus 
the period calculated from the mean energy, (E), which has also been eval- 
uated numerically. From this figure it appears that our qualitative argument 
of a trajectory stabilized around the Hamiltonian orbit, corresponding to the 
average energy, is well supported by the numerical simulations. 

In order to proceed further, we consider the probability distribution for the 
amplitude variables, P(Ai, A 2 , A 3 ; t) which obeys a Fokker-Planck equation 
[19]. For a general dynamics of the type given by Eq. (8), it is possible to 
show [15] that the stationary probability distribution for the Aj variables is 
given by 

P st (A u A 2 , A 3 ) = Z- 1 exp[-V(A u A 2 , A 3 )/e\, (50) 
Z = J dA 1 dA* 1 dA 2 dA* 2 dA 3 dA* 3 e~ v/e 

whenever two conditions are satisfied: 

a) Orthogonality condition (12). 

b) The residual dynamics [nonrelaxational part of (8)] is divergence free: 

Ova , . 

VMr°- (51) 

In our case of the Busse-Heikes equations the orthogonality condition is sat- 
isfied for /j, — 0, 8 > 0, and (51) is satisfied independently of /x and 5. For 
/i > this is no longer true but we expect that for small ji a relation similar 
to (50) would be valid if we replace V by a function $ that differs from V in 
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terms that vanish for vanishing /i. Using this probability distribution, one can 
compute the average value of the variable E as: 



(E) = Z- 1 J 



dA l dA* 1 dA 2 dA* 2 dA 3 dA* 3 £exp[-$/e]. 



(52) 



We take the crude approximation $ = V and, after a change of variables to 
amplitude and phase, the mean value of the energy can then be computed as: 



where V and E are given in terms of the variables R\, R 2 , R3 in Eqs. (9) and 
(27), respectively. In the case /i — (for which the above expression is exact) 
we obtain the value (E) = 1/60, independent of e, and T = T((E) = 1/60) m 
6.4467/5. 

In the case /i > 0, the above integral can be performed by means of a steepest 
descent calculation, valid in the limit e — » 0, where it yields the asymptotic 
behavior (E) — > (e/fi) 2 . The mean period can now be computed, in this limit 
of small e, using (43), with the result that the period, as a function of the 
system parameters 5, /i, e, behaves as: 



a relation that is expected to hold in the limit of small e and for small values 
of fi. The dependence with e is the same than the one holding for the mean 
first passage time in the decay from an unstable state [15] and also follows 
from the general arguments of [8]. In Fig. 8 we show that there is indeed a 
linear relation between the period computed in the numerical simulations and 
5 _1 ln(///e), as predicted by the above formula, although the exact prefactor 
3 is not reproduced. We find it remarkable that, in view of the simplifications 
involved in our treatment, this linear relation holds for a large range of values 
for the parameters /i, 5 and e. 



4 Spatial-dependent terms 

Tu and Cross [5] have proposed an alternative explanation for the stabilization 
of the period without the necessity of the inclusion of the noise terms: they 
modify the Busse-Heikes equations by considering two-dimensional amplitude 
fields, Aj(r,t), and including terms accounting for the spatial variation of 
those fields: 



(E) 



Jo 00 dR, Jo 00 dR 2 Jo 00 dR 3 Eexp[-V/e] 
Jo 00 dR l Jo 00 dR 2 / °° dR 3 expf-y/e] 



(53) 



T(e,/i,<5)^ln(/i/e), 



(54) 
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d t A 1 = C 1 A 1 + A 1 [l- 
d t A 2 = £ 2 A 2 + A 2 [l- 
d t A 3 = C 3 A 3 + A 3 [l- 



A 1 \ 2 -(l + f x + 5)\A 2 
A 2 \ 2 -(l + n + 5)\A 3 
A 3 \ 2 - (1 + ^ + 5)1^ 



(l + f i-6)\A 3 \% 

(l + ^-5)|A| 2 ], (55) 

(l + »-6)\A 2 \ 2 }. 



Here Cj (j — 1, 2, 3) are linear differential operators. Two main classes of op- 
erators can be considered: isotropic and anisotropic. Whereas a multiple scale 
analysis of the convective instability usually leads to anisotropic terms, the 
isotropic terms are often justified for the sake of mathematical and numer- 
ical simplicity. There are also the natural choice in problems of population 
dynamics [10]. The simplest isotropic terms are the Laplacian operators: 

4 = V 2 , J = 1,2,3. (56) 



Two types of anisotropic terms have been proposed for similar fluid problems 
in the literature: (i) the Newell- Whitehead-Segel (NWS) terms [20] and (ii) 
the Gunaratne-Ouyang-Swinney (GOS) terms [21,22]. Without altering the 
essentials of the problem, both NWS and GOS terms can be further simplified 
leading to second-order directional derivatives along three directions with a 
relative orientation of 60° [5,23]: 

Cf = (e r V) 2 , J = 1,2,3, (57) 



which are the only anisotropic terms considered henceforth. These are more 
tractable numerically and will be used to compute the alternating period as 
explained below. 

In this section we will compare the dynamical evolution corresponding to each 
one of the isotropic and anisotropic spatial dependent terms presented before, 
Eqs. (56) and (57), respectively. 

Common to all of them is that, as in section 2, we can recast system (55) into 
the form: 



a t A j (r,t) = -^ + 5v j , J = 1,2,3, 



6A* 



(58) 



where JF BH is a real functional of the fields given by: 



Fbu[A 1 ,A 2 ,A 3 ]= J rfr^Ql^A,! 2 - Kf + i|A/) + 
{1 + ^A^A^ 2 + \A 2 \ 2 \A 3 \ 2 + \A 3 \ 2 \A,\ 



(59) 
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and the functions Vj are given by (10). 



As in the zero-dimensional case of sections 2 and 3, 5 = entails a relax- 
ational gradient type dynamics and JF BH acts as a Lyapunov functional that 
decreases monotonically with time. Since this potential is minimized by homo- 
geneous solutions (because the spatial-dependent term gives always a positive 
contribution) the stationary solutions (and their stability) in the case 5 = 
are the same as in the zero-dimensional case. Unfortunately, the orthogonality 
condition 

*E/ dr ^ w ; + c - c - = ° (so) 



is not trivially satisfied in the case /x = for any of the spatial dependent 
terms mentioned before, and the dynamical equations can not be reduced for 
\i = as in the zero-dimensional case. 

In general, for < 5 < //, when the amplitudes grow from random initial 
conditions around Aj = 0, j = 1, 2, 3, we expect the formation of interfaces 
between the roll homogeneous states. Those interfaces move due to curvature 
and non-potential (5 > 0) effects. Moreover, the fact of dealing with three 
fields allows the formation of vertices, or points at which the three amplitudes 
take the same value. In the potential case, 5 = 0, the interface motion is 
such that a final state in which a unique roll solution fills the whole space is 
obtained (a process defined as "coarsening"). On the other hand, the nonpo- 
tential dynamics induces the rotation of front lines around vertices giving rise 
to the formation of rotating spiral structures [24]. Similar structures have been 
observed in other three competing species systems, such as lattice voter mod- 
els [25]. For small values of ji, the interfaces are wide (it can be shown that an 
interface varies over a length scale of order l/y/Ji) and the density of vertices is 
low. For large \i the interfaces are sharp and the density of vertices increases. 
The exact shape of the spirals depends upon the spatial-derivative terms used. 
With the isotropic terms, Eq. (56), interface propagation follows the normal 
direction at each point so that closed domains have spherical shape and spiral 
structures are close to Archimedes' spirals. On the other hand, for anisotropic 
spatial derivatives, Eq. (57), interface propagation no longer follows the normal 
direction and closed domains stretch or collapse along preferential directions 
so that they adopt an elliptic shape rather than a spherical one. 

An important effect is that the rotation of interfaces around vertices, driven 
by nonpotential effects, prevents the system from reaching a single roll solu- 
tion filling the whole space, even outside the Kuppers-Lortz instability region, 
i.e. for 5 < fi [26]. While this is true both for isotropic and anisotropic deriva- 
tives, the dynamical mechanism that prevents this coarsening is different for 
isotropic and anisotropic terms. For the isotropic terms, vertices of opposite 
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sense of rotation annihilate initially with each other if located closer than a 
critical distance d c ~ After a transient time in which vertices are formed, 
they place each other outside the range of effective attraction of other vertices 
so that their number is essentially constant, thus preventing coarsening. For 
the anisotropic terms, two interfaces associated with the same vertex (and 
thus rotating in the same sense) may collide and generate continuously new 
vertices which, in turn, annihilate against each other again preventing coars- 
ening outside the instability region. A consequence of interface motion is that 
a fixed point in space sees a change of the dominating amplitude. This alterna- 
tion change is essentially periodic in time and presents a characteristic period 
which has nothing in common with the Kiippers-Lortz instability mechanism 
in the bulk. Therefore, the period associated to this rotation is continuous at 
5 — //, the instability point. 

Before discussing what happens when this interface motion appears together 
with the instability in the bulk, we mention that for the isotropic terms it 
is possible to establish an analytical result concerning the front and spiral 
motion. In this case, using the fact that interfaces move in the normal direction 
to each point, it is possible to show that the rotation angular velocity of the 
interfaces around an isolated vertex scales, for small 5, as uo ~ 5 2 [26]. This 
predicts that, for the isotropic derivatives, the average period in a fixed point 
of space coming from the rotating spirals scales as (T) ~ 5~ 2 . 

As mentioned above, the mechanism of front motion due to the nonpotential 
effects coexists with the Kiippers-Lortz bulk instability. We will show in the 
remaining of the section some results that follow, mainly, from a numerical in- 
tegration of Eqs. (55) in two spatial dimensions. It appears from the numerical 
simulations that the behavior beyond the instability point (for 5 > fi) depends 
strongly on the type of spatial derivatives used as well on the magnitude of 
the parameter \i. We discuss first each type of derivatives separately. 

Isotropic derivatives: For // small, the bulk instability is such that the 
intrinsic Kiippers-Lortz period stabilizes to a statistically constant value. In a 
given point of space, we can see that the dominant amplitude changes due both 
to invasion from a rotating interface and a new amplitude growing inside the 
bulk. We give evidence of this combined mechanism in Fig. 9 where we have 
used the value /i — 0.1 and we present representative configurations inside and 
outside the instability region. 

For higher values of /i, the Kiippers-Lortz intrinsic period in the bulk is ob- 
served to increase with time. This is the same phenomenon that occurs in 
the zero-dimensional model without noise, see section II. Therefore at long 
times the Kiippers-Lortz period is so large that we only see rotating interfaces 
around vertices, just like below the instability point. The two images of the 
upper row in Fig. 10 show domain configurations at long times for ji = 2.5, 
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below (5 = 2) and beyond (S = 3.5) the Kiippers-Lortz instability point in 
the case of the isotropic terms. Apart from the typical size of the domains, it 
appears that there is no qualitative difference between them. The period of 
alternating amplitudes is entirely dominated by front motion. 

Anisotropic derivatives: Both for small and large /x, in the Kiippers-Lortz 
regime, 5 > /x, we observe, in addition to the front motion, domains of one 
phase emerging in the bulk of other domains; this is seen at all times, indicating 
that, at variance with the isotropic derivative case, the period associated with 
the Kiippers-Lortz instability does not diverge with time. Evidence is given in 
Fig. 9 for /x = 0.1 and Fig. 10 for /x = 2.5, both figures showing results inside 
and outside the instability region. 

For small /x, in summary, the morphology of domains inside and outside the 
instability region turns out to be similar with both kinds of spatial depen- 
dent terms, Fig. 9. The alternating period for 8 > /x is dominated by the 
Kiippers-Lortz instability and is similar with isotropic and anisotropic spatial 
derivatives. This shows up in the fact that the period computed in a single 
point of space does not depend essentially of the type of derivatives used, as 
shown in Fig. 11a. 

For large /x, on the other hand, the morphology is different for isotropic and 
anisotropic terms. For the isotropic ones, spiral rotation dominates the dy- 
namics because of the very large period associated with the bulk instability. 
For the anisotropic terms, both front motion and bulk instability are present. 
Finally in Fig. lib (large /x) we show how the alternating period changes when 
going through the Kiippers-Lortz instability. We first note that the period does 
not vanish in the stable regime (S < /x). In this regime it is entirely due to 
front and spiral motion. For isotropic derivatives the period changes smoothly 
through the point 5 — /x. This supports the fact that the period is still given 
by front motion for 8 > /x. On the contrary, for anisotropic derivatives a jump 
in T is observed at 5 — /x. In the Kiippers-Lortz unstable regime and for 
anisotropic derivatives, T is determined by a combination of bulk instability 
and front motion. 



5 Conclusions 

We have analyzed the Busse-Heikes equations for Rayleigh-Benard convection 
in a rotating fluid. For the situation of spatial-independent amplitudes, a case 
previously analyzed by May and Leonard, we find a Lyapunov potential that 
allows us, for /x = 0, to split the dynamics into a relaxational plus a residual 
part. Since the residual dynamics is Hamiltonian, we are able to give explicit 
relations for the time variation of the amplitudes and to compute the period 
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of the orbits as a function of the energy, which, in turn, is a function of 
initial conditions. For /i > we extend the previous picture by using an 
adiabatic approximation in which the energy slowly decreases with time. This 
allows us to compute the variation of the alternation period between the three 
modes in the Kuppers-Lortz instability regime. We next consider the effect 
of fluctuations and show how noise can stabilize the mean period to a finite 
value. By using the Lyapunov potential employed in the deterministic case, we 
can deduce an approximate expression that yields the period as a function of 
the system parameters, //, 5 as well as a function of the noise intensity e. The 
conclusion is that the period increases logarithmically with decreasing noise 
intensity, a result that is well confirmed by numerical simulations 

The two-dimensional version of this problem exhibits rather different dynami- 
cal behavior grossly dominated by vertices where three domain walls meet and 
which have no parallel in lower dimensional systems. The rotation of interfaces 
around vertices is driven by nonpotential effects and this inhibits coarsening 
for sufficiently large systems. We investigated the influence on the dynamics of 
the type of spatial dependent terms. For small values of the parameter /i, the 
morphology of domains inside the Kuppers-Lortz region turns out to be similar 
for both isotropic and anisotropic spatial derivatives. The alternating period 
is dominated by the Kuppers-Lortz instability and is similar for both kinds 
of spatial-dependent terms. For large /i, on the contrary, the morphology of 
patterns as well as the alternating mean period are different for isotropic and 
anisotropic terms. While the intrinsic period of the instability diverges with 
time with isotropic derivatives, it saturates to a finite value in the anisotropic 
case. 

We acknowledge financial support from DGESIC (Spain) projects numbers 
PB94-1167 and PB97-0141-C02-01. 
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Fig. 1. Stability regions for the Busse-Heikes dynamical system (2). The region 'H' 
is where the hexagon solution (three equal amplitudes) is stable. In the 'R' region, 
the three roll solutions are stable, and in region 'KL' there are no stable fixed points. 




Fig. 2. Dynamics for ji = in the variables Ri, R2, R3 for two different initial con- 
ditions. After a transient time of order 1 the motion is on the plane R1+R2 + R3 = 1 
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Fig. 3. Time evolution of amplitudes in the case 5 = 1.3, fi = 0. After a transient 
time of order 1, the three variables R±, R2, R3 vary periodically in time. The lines 
are the theoretical predictions that come from Eqs. (40,41,42). 
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Fig. 4. Time evolution of amplitudes in the case 8 = 1.3, /i = 0.1. The characteristic 
alternation time of the three variables R\, R2, R3 increases with time. Notice that 
the envelope of the amplitudes approaches one asymptotically, and that their sum, 
R\ + R2 + Rz is approximately equal to 1. 
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Fig. 5. Time evolution of the period, defined as the time it takes a given amplitude 
to cross the reference level Rj = 0.5 plotted versus time for several values of 5 and 
(x. We also plot straight lines with slopes -j- as predicted by Eq. (46). 




Fig. 6. Time evolution of amplitudes in the presence of noise for 5 = 1.30, [i = 0.1, 
e = 10~ 7 . In this case, the motion is such that the time interval between dominations 
of a single mode fluctuates around a mean value (compare with the equivalent 
deterministic case shown in Fig. 4). 
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Fig. 7. Plot of the average period (T) plotted versus the theoretical value T((E)) 
computed using the value of (E) obtained numerically. For each value of fj, and S 
(same symbols meaning than in Fig. 5) we use values of e ranging from e = 10~ 2 to 
e = 10~ 7 . 




Fig. 8. Average period, (T), plotted as a function of <5 _1 log( / u/e) in order to check 
the predicted linear relation (54). The straight line is the best fit and has a slope of 
1.73. Same symbols meanings than in Fig. 5 and values of e ranging from e = 10~ 2 
to e = 10~ 7 . 
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Fig. 9. Four snapshots at long times corresponding to a numerical simulation 
of the Busse-Heikes model [eq. (55)] with isotropic (£j = V 2 ) and anisotropic 
{C"j = (ej • V) 2 ) spatial derivatives. Parameter values are: fi = 0.1 and 5 = 0.05 
(1.3) outside (inside) the Kiippers-Lortz instability region. 
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Fig. 10. Four snapshots at long times corresponding to a numerical simulation 
of the Busse-Heikes model [eq. (55)] with isotropic (£j = V 2 ) and anisotropic 
{Cf = (etj ■ V) 2 ) spatial derivatives. Parameter values are: jjl = 2.5 and 5 = 2 (3.5) 
outside (inside) the Kiippers-Lortz instability region. 



(a) |LL = 0.1 



(b) |a = 2.5 
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Fig. 11. Inverse of the alternating mean period as a function of 5 2 for 
the two-dimensional Busse-Heikes model with isotropic and anisotropic spa- 
tial-dependent terms. We have chosen the coordinates in order to emphasize the 
linear relation between the inverse of the period and 5 2 valid for small S [26]. Each 
plot corresponds to a different value of the parameter /i. The Kiippers-Lortz insta- 
bility takes place at the right of the vertical dotted lines. 
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